The Journal of Chemical Physics

Molecular dynamics (MD) simulations are a powerful tool for investigating biomolecular dynamics underlying biological functions. However, the accessible spatiotemporal scales of conventional all-atom simulations remain limited by high computational costs. Coarse-graining reduces these costs by decreasing the number of interaction sites and enabling longer timesteps. In extreme cases, proteins are represented as single spherical particles; while such approximations facilitate cellular-scale simulations, they often sacrifice essential structural information, such as molecular shape and interaction anisotropy. Here, we present CGRig, a rigid-body protein model with residue-level interaction sites designed for long-time, large-scale simulations. In CGRig, each protein is treated as a single rigid-body embedding residue-level interaction sites. Its translational and rotational motions are described by the overdamped Langevin equation incorporating a shape-dependent friction matrix. Intermolecular interactions are calculated using G[o]-like native contact potentials, Debye-Huckel electrostatics, and volume exclusion. We validated that CGRig accurately reproduces the translational and rotational diffusion coefficients expected from the friction matrix for an isolated protein. For dimeric systems, the model successfully maintained native complex structures. Furthermore, two initially separated proteins converged into the correct complex with an association rate consistent with all-atom simulations. Notably, CGRig achieved a simulation performance exceeding 17 s/day for a 1,024-molecule system. These results demonstrate that CGRig provides an efficient framework for simulating protein assembly while retaining residue-level interaction specificity, making it a valuable tool for investigating large-scale biomolecular self-assembly.

There is a growing demand for molecular dynamics simulations to explore longer timescale behavior of giant protein-DNA complexes such as chromatin. To address this need, we extended OpenCafeMol, a GPU-accelerated residue-level coarse-grained molecular dynamics simulator originally developed for proteins and lipids, to support 3SPN.2 and 3SPN.2C DNA models. We also implemented a hydrogen-bond-type many-body potential to model DNA-protein interactions more accurately. To further improve computational efficiency, we introduced a localized scheme for calculating base-pairing and cross-stacking interactions. Benchmark tests show that OpenCafeMol on a single GPU achieves up to 200-fold speed-up for DNA-only systems and up to 100-fold speed-up for DNA-protein complexes compared to CPU-based simulations. To demonstrate the capability of our implementation for long-timescale biological processes, we simulated an archaeal SMC-ScpA complex undergoing DNA translocation via segment capture (a proposed mechanism for DNA loop extrusion) in the presence of a DNA-bound obstacle. We observed continuous captured-loop growth accompanied by obstacle bypass within the segment capture framework.

Accurate ion force fields are essential for molecular dynamics simulations of biomolecular systems, particularly in combination with modern water models such as OPC. While OPC water improves the description of bulk water and biomolecules, the transferability of existing ion force fields to this model remains an open question. Here, we systematically assess the transferability of monovalent and divalent ion force field parameters (Li+, Na+, K+, Cs+, Mg2+,Ca2+, Sr2+, Ba2+, Cl- and Br-) to OPC water by comparing single-ion and ion-pairing properties with experimental data. Our analysis reveals that no single literature parameter set provides accurate results for all ions when directly transferred to OPC water. We hence introduce the MS/G-LB(OPC) force field, which combines Mamatkulov-Schwierz-Grotz cation parameters with Loche-Bonthuis anion parameters. MS/G-LB(OPC) reproduces hydration free energies, first-shell structural properties and activity derivatives at low salt concentrations. Our results demonstrate that transferring ion parameters to OPC can lead to significant and ion-specific deviations from experimental data, making careful validation essential. At the same time, the systematic transfer and combination of ion parameters from existing force fields can provide a practical and computationally efficient alternative to full reparameterization. MS/G-LB(OPC) is available at https://git.rz.uni-augsburg.de/cbio-gitpub/opc-ion-force-fields.

Biomolecular condensates formed by complex coacervation of highly charged proteins provide a powerful framework to understand how microscopic interactions give rise to macroscopic material properties. Atomistic molecular dynamics simulations provide detailed insights but remain limited in accesing the spatio-temporal scales relevant for condensate behavior. Here, we use the residue-level coarse-grained Mpipi-Recharged model to investigate condensates formed by ProT and positively charged partners, including histone H1, protamine, poly-lysine, and poly-arginine. Material properties, in this context, provide a stringent experimental benchamark for coarse-grained models. Our model reproduces salt-dependent phase behavior, protein binding affinities, and sequence-specific stability trends in agreement with in vitro experiments, despite the fact that material properties were not included in the model parametrization. We then establish a direct link between protein dynamics and macroscopic material properties by quantifying monomeric diffusion, conformational reconfiguration, and translational mobility within the dense phase, and relating these to condensate viscosity. By comparing dynamics across dense and dilute phases, we uncover a pronounced length scale-dependent behavior. While residue-level binding and unbinding events remain equally fast in both phases, protein reconfiguration time and self-diffusion are significantly slowed down within the condensates. This decoupling reveals how fast intermolecular interactions coexist with slow mesoscale condensate dynamics depending on the molecular length scale. Together, our results establish a predictive framework that links encoded sequence intermolecular forces and multiscale dynamics to the emergent material properties of complex biomolecular condensates.

The interactions of droplets and filaments can lead to mutual deformations and complex combined behavior. Such interactions also occur within the cell, where biomolecular condensates, distinct liquid phases often composed of proteins, have been observed to structure and affect the organization of the cytoskeleton. In particular, biomolecular condensates have been shown to undergo characteristic deformations when cytoskeletal filaments are fully embedded within them. However, a full understanding of the underlying physical mechanisms is still missing. Here, we combine experiments with coarse-grained molecular dynamics simulations and analytical models to uncover the physical mechanisms that define emerging shapes of droplets containing filaments. We find that the surface tension of the liquid phase and the bending energy of the filament(s) suffice to accurately capture emerging shapes if the length of the filament is small compared to the liquid volume. As the volume fraction of filament(s) increases, wetting effects become increasingly important, setting physical constraints within which surface and bending energies compete to define the droplet shapes. We find that mutual deformations of condensate and filament extend accessible shapes beyond classical stability considerations, leading to structuring and entrapment of contained filaments. Shape deformations may further affect ripening dynamics that favor certain geometries. Our findings provide a physical framework for a better understanding of the possible roles of biomolecular condensates in cytoskeletal organization.

We develop a self-consistent free-energy framework in which membrane shape and osmotic pressure are determined simultaneously in a finite reservoir by minimizing bending elasticity and solute entropy. Solute conservation makes osmotic pressure a thermodynamic variable rather than an externally prescribed parameter, producing a nonlinear coupling between membrane mechanics and solvent entropy. This coupling modifies the classical stability condition for spherical vesicles: instability emerges from global free-energy competition rather than the linear Helfrich stability criterion. The resulting critical pressures differ by orders of magnitude from Helfrich predictions and agree with simulations for small and large unilamellar vesicles. The framework is relevant to cellular environments involving biomolecular condensate confinement as well as synthetic vesicles and the development of osmotic-pressure-driven encapsulation platforms.

We present a new formulation of the low-field effect (LFE) in spin-correlated radical pairs based on a zero-field singlet-triplet basis for the isotropic spin Hamiltonian. The aim is to provide a description that is both formally rigorous and mechanistically transparent, especially in the regime of weak magnetic fields such as the geomagnetic field. For the standard model radical pair containing a single spin [Formula] nucleus, we show that the conventional singlet-triplet basis obscures the distinct dynamical roles of the hyperfine and Zeeman interactions. In the zero-field S-T basis, by contrast, the mechanism separates cleanly: isotropic hyperfine coupling mixes singlet-doublet and triplet-doublet states, whereas the weak-field Zeeman interaction mixes triplet-quartet and triplet-doublet states without directly introducing an additional singlet-triplet coupling. The LFE is therefore revealed as a sequential process in which a weak field unlocks access from a triplet-only manifold to a singlet-accessible triplet manifold, from which hyperfine-driven singlet-triplet interconversion can occur. We then generalize this picture to radical pairs with arbitrary isotropic hyperfine structures by identifying maximal, interior, and, when present, minimal triplet-only manifolds in the zero-field spectrum. Finally, we introduce a practical blockwise dark-state recruitment measure for the triplet-only zero-field state space made singlet-accessible by a weak field, and show how this quantity depends on hyperfine symmetry, including the effects of equivalent nuclei. The resulting framework provides both a simple physical picture of the LFE and a general route to estimating its structural upper bound for arbitrary radical pairs.

Phase separated condensates are recognized as a ubiquitous mechanism of spatial organization in cell biology. Biophysical modeling of condensates provides critical insights into the dynamics and functions of these subcellular structures that are difficult to extract via experiments. Here we present an efficient computational pipeline, CASPULE (Condensate Analysis of Sticker Spacer Polymers Using the LAMMPS Engine), to simulate and analyze the biological condensates made of sticker-spacer polymers. CASPULE implements a unique force field that combines traditional Langevin dynamics with a "detailed balance proof" protocol for single-valent bond formation between stickers. This framework allows us to study the non-trivial biophysics that emerge out of the single-valent sticker interactions coupled with the effect of separation in energetic contribution by stickers and spacers. We provide detailed documentation on how to setup the simulation environment, perform simulations and analyze the results. Through case studies, we highlight the utility and efficacy of our pipeline. Importantly, we provide statistical parameters to characterize the cluster size distribution often observed in biological systems. We envision this tool to be broadly useful in decoding the interplay of kinetics and thermodynamics underlying the formation and function of biological condensates.

The apparent pKa of ionizable lipids in lipid nanoparticles (LNPs) is a key determinant of RNA encapsulation during formulation and endosomal release after cellular uptake. However, it is difficult to predict the effective pKa of a given ionizable lipid solely from its solution pKa, because it is sensitive to the membranes composition, as well as solution conditions such as the salt concentration. We developed a simple continuum electrostatics model, based on Gouy-Chapman theory, to predict the shift in effective pKa for ionizable lipids in lipid bilayers as a function of salt concentration and membrane composition. We derive equations for the surface potential and fraction of lipids charged, which are solved self-consistently as a function of solution pH to extract the titration curve and effective pKa. The model shows that the shift in effective pKa is largest when the concentration of titratable lipid is high, and the effect is diminished by increasing salt concentration. We provide a python implementation of the model and an interactive notebook that will allow users to further easily explore the predicted pKa shifts as a function of formulation variables.

Cerebrospinal fluid (CSF) is a Newtonian fluid that bathes the brain and spinal cord and oscillates in response to the physiological periodic changes in brain volume, of which the cardiac cycle is a major driver. Understanding this motion is essential for clarifying its contribution to solute transport, waste clearance, and drug delivery. In this work, we study oscillatory and steady streaming flow in the cranial subarachnoid space using a lubrication-based theoretical framework. The model represents the cranial CSF compartment as a thin fluid layer bounded internally by the brain surface and externally by the dura, driven by time-dependent brain surface displacements. We first derive simplified governing equations for flow over an arbitrary smooth sphere-like brain surface and obtain analytical solutions for an idealised spherical geometry with uniform displacements. We then incorporate realistic displacement fields reconstructed from MRI measurements in healthy subjects and solve the reduced equations numerically. The results show that oscillatory forcing produces a steady streaming component that may enhance solute transport compared with diffusion alone. This work provides a mechanistic description of the flow generated by physiological brain motion and highlights the potential presence of steady streaming in cranial subarachnoid fluid dynamics.

Protein kinases regulate signaling by recognizing short sequence motifs, and how these motifs bind influences both specificity and therapeutic strategies that target kinase pathways. Peptide-based inhibitors that engage substrate-recognition regions are attracting interest, but designing them requires an understanding of how a flexible peptide approaches and settles into the bound pose. Traditional studies have focused on the bound pose and affinities, whereas the steps that link the initial encounter with the bound pose have been explored less thoroughly because the relevant intermediates are too short-lived to capture experimentally and evolve on timescales that standard molecular dynamics cannot readily access. Here, we focused on Abl kinase and Abltide, the experimentally identified optimal substrate peptide for Abl kinase, and examined the sequence of events linking initial encounter to the bound pose using two-dimensional replica exchange (gREST/REUS), which selectively enhances flexibility in the peptide and its binding interface while also sampling progression along a distance coordinate. The resulting simulations yielded a detailed binding landscape, revealing five distinct encounter regions outside the substrate-binding site and six intermediate states that may connect the initial approach to the bound pose. Some encounter regions and intermediate states participate in the dominant binding pathways. During this process, EF/G/{beta}11 hydrophobic patch, together with G helix negative patch, plays a central role in guiding Abltide toward the substrate-binding site. These findings provide mechanistic insight into substrate recognition by protein kinases and offer a foundation for the rational design of peptide-based inhibitors.

Understanding how cells migrate through confined environments is crucial for elucidating fundamental biological processes, including cancer invasion, immune surveillance, and tissue morphogenesis. The nucleus, as the largest and stiffest cellular organelle, often limits cellular deformability, making it a key factor in migration through narrow pores or highly constrained spaces. In this work, we introduce a geometric surface partial differential equation (GS-PDE) model in which the cell plasma membrane and nuclear envelope are described as evolving energetic closed surfaces governed by force-balance equations. We replicate the results of a biophysical experiment, where a microfluidic device is used to impose compressive stresses on cells by driving them through narrow microchannels under a controlled pressure gradient. The model is validated by reproducing cell entry into the microchannels. A parametric sensitivity analysis highlights the dominant influence of specific parameters, whose accurate estimation is essential for faithfully capturing the experimental setup. We found that surface tension and confinement geometry emerge as key determinants of translocation efficiency. Although tailored to this specific setup for validation purposes, the framework is sufficiently general to be applied to a broad range of cell mechanics scenarios, providing a robust and flexible tool for investigating the interplay between cell mechanics and confinement. It also offers a solid foundation for future extensions integrating more complex biochemical processes such as active confined migration.

Recently, a number of tools have been released that generate ensembles of protein structures based on artificial intelligence (AI) approaches. Although ensembles generated by the tools differ significantly, we demonstrate a computational path to harmonizing the various outputs under a stationary condition using two complementary physics-based approaches. In the first stage, the AI ensemble is used to seed a weighted ensemble (WE) simulation, promoting relaxation toward the steady state. In the second stage, trajectory segments generated by WE are reweighted to steady state using the recently developed RiteWeight (RW) algorithm. We applied this approach to generate an atomically-detailed equilibrium ensemble of unliganded adenylate kinase conformations, starting from ensembles produced by three AI tools: AFSample2, ESMFlow-PDB (trained from PDB structures), and ESMFlow-MD (trained from molecular dynamics simulation data). Dramatic differences in the AI-generated ensembles are largely erased during the WE-RW process, yielding a consistent description of the equilibrium ensemble for a given force field.

Reactive oxygen species generated during inflammation can oxidize viral envelope lipids, with outcomes ranging from modulated infectivity to viral inactivation. For SARS-CoV-2, the molecular mechanisms by which membrane lipid oxidation influences spike protein anchoring remain poorly understood. We use all-atom molecular dynamics (MD) simulations to quantify how graded oxidation of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) affects the anchoring of the SARS-CoV-2 spike transmembrane (TM) region in an endoplasmic-reticulum-Golgi intermediate compartment (ERGIC)-like multicomponent membrane. Viral envelopes containing 0, 25, 50, 75, and 100% oxidized POPC (PoxnoPC) corresponding to 0 - 55% oxidation of all PO-type phospholipids were simulated with the spike TM helix and cytoplasmic tail embedded in a POPC/POPE/POPI/POPS/cholesterol mixture. Steered MD and umbrella sampling were used to calculate the potential of mean force (PMF) for extracting the TM+CT region along the membrane normal. Partial oxidation (25 - 75% POPC) produced reductions in the detachment barrier that were not statistically distinguishable from the native system within the sampling uncertainty, whereas full POPC oxidation lowered the anchoring free energy by about 23% (from 606 {+/-} 39 to 464 {+/-} 38 kJ mol-1), indicating that oxidation of roughly half of the glycerophospholipids can measurably weaken spike-membrane coupling. Despite this reduction, the remaining barrier (about 180kBT ) is still large, suggesting that oxidation alone may be insufficient for spontaneous spike detachment and likely acts synergistically with mechanical forces during fusion or immune engagement. Analysis of acyl-chain order parameters, area per lipid, membrane thickness, number-density profiles, and lateral lipid clustering reveals that POPC peroxidation decreases lipid order, thins and softens the bilayer, and disrupts cholesterol-stabilized clusters that refer to large cooperative lipid assemblies (>10 lipids) identified via RDF-based clustering. These oxidation-induced changes reduce hydrophobic matching around the TM helix and facilitate its extraction from the viral envelope. Our results provide a mechanistic link between lipid peroxidation, membrane nanostructure, and spike anchoring, supporting lipid oxidation for example during cold atmospheric plasma or ozone treatment as a physically grounded contributing antiviral mechanism against SARS-CoV-2.

Cytokine-mediated communication is a central mechanism by which immune cells coordinate activation, differentiation and proliferation. While mechanistic reaction-diffusion models provide detailed descriptions of cytokine secretion and uptake at the cellular scale, their computational cost limits their applicability to large and densely packed cell populations. Previously employed approximations of cytokine diffusion fields rely on assumptions that neglect the influence of cellular geometry and volume exclusion. In this work, we study a macroscopic description of cytokine diffusion and reaction dynamics based on homogenization techniques, rigorously linking microscopic reaction-diffusion formulations to effective continuum models. The resulting homogenized equations replace discrete responder cells with a continuous density, while retaining essential features of cellular uptake and excluded-volume effects. Further, we show that in regimes with approximate radial symmetry, classical Yukawa-type solutions emerge as limiting cases of the homogenized model, provided appropriate correction factors are included. Overall, our approach allows efficient multiscale modeling of cytokine signaling in complex immune-cell environments.

It is an essential element of mechanobiology to measure the forces of biological cells. In microparticle traction force microscopy, they are inferred from the deformation of elastic microparticles. Two complementary variants have been introduced before: the volume method, which reconstructs surface stresses from the displacements of fiducial markers embedded inside the particles, and the surface method, which infers stresses directly from the deformation of the particle surface. However, a systematic comparison of the two methods has been lacking. Here, we quantitatively compare both approaches using simulated traction fields representing biologically relevant loading scenarios. We find that the surface method consistently reconstructs traction profiles with substantially lower errors than the volume method, which suffers from displacement tracking and stress calculation at the surface. At high noise levels, however, the performance gap becomes smaller. To compare the performance of the two methods in a realistic experimental setting, we developed DNA-based hydrogel microparticles equipped with both fluorescent surface labels and embedded fluorescent nanoparticles, enabling the direct comparison of the two methods within the same system. Compression experiments produced traction profiles consistent with Hertzian contact mechanics and confirmed the trends observed in the simulations. While our computational workflow establishes a framework to apply both methods, our experimental workflow establishes DNA microparticles as versatile and biocompatible probes for measuring cellular forces.

The rational design of biomolecule immobilization strategies requires molecular-level understanding of how surface properties, tethering geometry, and structural dynamics jointly influence stability and function. Recently, coarse-grained molecular dynamics simulations based on the Martini force field have emerged as an efficient framework for studying enzyme-surface interactions. However, the reproducible construction of immobilized systems with controlled orientations remains technically challenging, usually involving multiple computational tools. Here we present MartiniSurf, an open-source command-line framework for the preparation of protein and DNA systems immobilized on solid supports within the Martini paradigm. MartiniSurf integrates automated structure retrieval and cleaning, coarse graining via tools from the Martini force field software ecosystem, customizable surface generation, and biomolecule orientation based on user-defined anchoring residues, producing complete GROMACS-ready simulation systems. The framework supports both implicit restraint-based anchoring and explicit linker-mediated immobilization, including surfaces functionalized with user-defined ligands or linker-like moieties, enabling representation of mono- and multivalent attachment geometries at different modeling resolutions. Structure-based G[o]Martini potentials can be incorporated for proteins, while DNA systems are modeled using Martini 2. Optional substrate insertion, pre-coarse-grained complex handling, and automated solvation and ionization further extend system flexibility. By integrating these components into a unified workflow, MartiniSurf enables systematic and high-throughput in silico exploration of surface-tethered biomolecules and provides a robust computational platform for rational immobilization studies. TOC Graphic O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=146 SRC="FIGDIR/small/714767v1_ufig1.gif" ALT="Figure 1"> View larger version (45K): org.highwire.dtl.DTLVardef@bc1ac4org.highwire.dtl.DTLVardef@1813b43org.highwire.dtl.DTLVardef@159b19borg.highwire.dtl.DTLVardef@19b60d6_HPS_FORMAT_FIGEXP M_FIG C_FIG

Understanding RNA structural dynamics is essential for elucidating its biological functions, and molecular dynamics (MD) simulations provide an important atomistic complement to experimental approaches. However, the predictive power of MD is fundamentally limited by the accuracy of the underlying empirical Force Fields (FFs), particularly in capturing the delicate balance of non-bonded interactions. Here, we present a systematic reparameterization strategy that replaces the external gHBfix19 hydrogen-bond (H-bond) correction potential with an equivalent set of NBfix Lennard-Jones modifications within a state-of-the-art RNA FF. Using a quantitatively converged temperature replica-exchange MD ensemble of the GAGA tetraloop, we employed a reweighting-based optimization protocol to derive NBfix parameters that reproduce the thermodynamic effects of the original gHBfix19 terms. Sequential optimization of individual gHBfix19 components proved essential to ensure stable and transferable parameter refinement. The resulting fully reformulated NBfix-based variant, termed OL3CP-NBfix19, was validated on a representative set of RNA motifs, including tetranucleotides, A-form duplexes, and tetraloops. Across all tested systems, its performance is comparable to that of the reference gHBfix19 FF. By embedding the H-bond corrections directly into the standard non-bonded framework, the NBfix formulation eliminates external biasing potentials, simplifies practical deployment, and reduces computational overhead. Beyond this specific reparameterization, our results demonstrate a practical workflow for translating targeted H-bond corrections into native FF terms for efficient biomolecular simulations.

Developing approaches to link structure and function is an ongoing challenge in computational and structural biology. Using a systems-level framework, we present here an analysis pipeline in a Python package, mdsa-tools, that constructs network representations of structures in a time series of trajectory frames from molecular dynamics (MD) simulations. Here, we demonstrate its use on a ribosomal subsystem. The subsystem is centered on the CAR interaction surface, a "brake pad" adjacent to the aminoacyl (A-site) decoding center that tunes protein translation rates. We leverage unsupervised learning algorithms to explore the conformational landscape of behaviors visited by two versions of the subsystem (brake-on and brake-off) that differ at the codon 3 adjacent to the A-site codon. Our network representations of MD frames embody H-bond interactions between all pairwise combinations of residues in the system. By utilizing per-frame vector representations of network edges, we can apply standard clustering and dimensionality reduction methods to explore behavioral differences between the brake-on and brake-off versions of the system. K-means clustering of frame vectors revealed a striking separation of the two system versions, consistent with principal components analysis (PCA) embeddings and Uniform Manifold Approximation and Projection (UMAP) embeddings. Dissection of K-means centroids and PCA loadings highlighted H-bond interactions between residue pairs in the ribosomes peptidyl site (P site), suggesting potential allosteric signaling across the subsystem. Author summaryWith the impressive development of computational algorithms to successfully simulate the dynamics of biological molecules over time, the exploration and incorporation of systems modes of analysis is a natural next step to begin to understand the molecular dynamics behaviors that emerge from these experiments. Following the approaches of classical molecular genetics, we used a "computational genetics" paradigm where we introduced changes (mutations) in potentially important residues, changing their identities or modifying their chemical properties, and asked how the dynamic system responded to these changes, viewing the simulations as a series of movie frames of the dynamic structure over time. Starting with network representations of each frames structure, where the nodes are residues, and the edges denote H-bond interactions between the residues, we used several unsupervised machine learning algorithms to uncover behavioral changes in the different mutated versions of the system. Applied to our ribosome neighborhood, this revealed unexpected changes in behavior at the ribosome peptidyl site (P site) in response to mutating mRNA residues on the other side of the aminoacyl site (A site) codon, suggesting long-range allosteric interactions across the neighborhood.

Alzheimers disease (AD) is characterized by the accumulation of amyloid-{beta} (A{beta}), yet the specific link between plaque burden and cognitive decline remains a subject of intense investigation. This paper presents a mathematical model that simulates the coupled dynamics of A{beta} monomers, soluble oligomers, and fibrillar species in the brain tissue. By modifying existing moment equations to include a dedicated conservation equation for A{beta} monomers, the model explores how various microscopic processes, such as primary nucleation, surface-catalyzed secondary nucleation, fibril elongation, and fragmentation, contribute to macroscopic disease progression. Central to this study is the concept of "accumulated neurotoxicity" as a surrogate marker of biological age, defined as the time-integrated concentration of soluble A{beta} oligomers. Unlike plaque burden, accumulated neurotoxicity cannot be reversed, and the harm it causes depends critically on the sequence of events that produced it. Numerical results demonstrate that while plaque burden and neurotoxicity both increase over time, their relationship is non-linear and highly sensitive to the efficiency of protein degradation machinery. Specifically, impaired degradation leads to a rapid advancement of biological age relative to calendar age. The model further identifies oligomer dissociation and fibril fragmentation as potential protective mechanisms that can counterintuitively reduce neurotoxic burden by diverting monomers away from the soluble oligomer pool. These findings provide a quantitative framework for understanding why individuals with similar plaque burdens may experience vastly different cognitive outcomes, underscoring the importance of targeting soluble oligomers early in therapeutic interventions.